Influence of the temperature on the depinning transition of driven interfaces 
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We study the dynamics of a driven interface in a two- 
dimensional random-field Ising model close to the depinning 
transition at small but finite temperatures T using Glauber 
dynamics. A square lattice is considered with an interface ini- 
tially in (ll)-direction. The drift velocity v is analyzed using 
finite size scaling at T = and additionally finite temperature 
scaling close to the depinning transition. In both cases a per- 
fect data collapse is obtained from which we deduce w 1/3 
for the exponent which determines the dependence of v on the 
driving field, v ~ 1 for the exponent of the correlation length 
and 8 « 5 for the exponent which determines the dependence 
of v on T. 

PACS: 68.35.Rh, 75.10.Hk, 75.40.Mg 



The motion of a driven interface in a random medium 
has attracted a large amount of interest recently because 
it is a challenging theoretical problem and because it oc- 
curs in many different areas of physics like fluid invasion 
of porous media [jlj , depinning of charge density waves ||] 
or field driven motion of a domain wall in a ferromagnet 
H , to name just a few. In these cases the disorder is time 
independent or quenched. For zero temperature there 
is a well defined critical driving force F c above which 
the interface moves while the interface gets trapped in 
metastable positions below it. This scenario has been 
observed in a variety of different works in the past Q. 

At finite temperature the depinning transition is 
smeared out and the interface can move also below F c . 
For very small driving forces a creep motion is expected 
which is governed by thermal activation while for driving 
forces close to the critical one a scaling behavior of the 
interface velocity v is expected to occur for small enough 
temperatures (see j| where the influence of the temper- 
ature on the velocity of sliding charge-density waves is 
discussed). 

The dynamics and the morphology of interfaces in sys- 
tems with quenched disorder have been investigated in 
the past within a variety of different models. Most stud- 
ies focus on equations of motion for the interface itself. 
Famous approaches are the Edwards- Wilkinson (EW) 
equation || or the Kardar-Parisi-Zhang (KPZ) equation 
H which have been studied extensively both with an- 
nealed and quenched disorder. A simple example for a 
disordered medium in which the motion of an interface 
can be studied is the random-field Ising model. In this 
case the interface is a domain wall separating regions 



of up and down spins. For this system it were Bru- 
insma and Aeppli H who argued with the assumption 
that the interface can be treated as an elastic membrane 
that its equation of motion at the depinning transition 
can be described by the EW equation with quenched dis- 
order. Their arguments are plausible but far from rig- 
orous. Therefore, it is also of great interest to tackle 
the full problem, i. e. to study the dynamics associated 
with the Hamiltonian of the random-field Ising model 
prepared initially with an interface. Under the influence 
of external forces the domain wall may begin to move in 
a way which depends on the strength of the disorder. A 
still unsolved problem is the scaling behavior of the ve- 
locity of the domain wall for small but finite temperature 
close to the critical driving field. This, of course, is an 
important issue also experimentally Q|. 

It is the purpose of this letter to elucidate this critical 
behavior of a domain wall in a two-dimensional random- 
field Ising model close to the pinning transition. Finite 
size scaling is used for the analysis of the zero tempera- 
ture behavior resulting in precise values of the exponents. 
For the first time we determine the critical exponent S as- 
sociated with the finite temperature smearing of the tran- 
sition which to the best of our knowledge is unknown at 
present. 

For a random-field Ising system in previous work || a 
cubic structure was considered with an interface between 
up and down spins initially parallel to one of the cubic 
axis of the lattice. An applied magnetic field H favor- 
ing the up spins energetically leaves the down spins in a 
metastable state. Under the influence of some dynamical 
rules the area occupied by metastable states shrinks, i. e. 
the interface starts to move. For the dynamics at zero 
temperature simple relaxation dynamics was assumed, 
i. e. a spin is flipped only if its energy is lowered. Addi- 
tionally it was assumed that only spins at the interface 
are allowed to flip in order to avoid spontaneous domain 
growth in the metastable state. The critical field was ap- 
proached from below, i. e. the driving field was increased 
in small steps until the pattern of flipped spins after an 
increase of the field spans the system. The field value for 
this to happen is the critical field H c . 

The geometry used in these investigations has the dis- 
advantage that even without impurities the interface is 
pinned for up to quite large fields. For space dimension 
d — 2, for instance, the driving field has to be larger than 
2J for the interface to move where J denotes the near- 
est neighbor exchange interaction. But if one spin of the 
interface is flipped all neighbors in the same line will flip 
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in the following updates. Thus for zero or very small dis- 
order the movement of the interface consists of a series 
of complete lines of spins which flip resulting in faceted 
growth. In this case pinning of the interface is not due to 
disorder but to the fact that any spin in the interface is 
locked by its four neighbors. Faceted growth can occur 
in any dimension for weak and bounded disorder but it 
depends on the structure of the lattice ||. 

By changing the geometry, i. e. considering a (11)- 
interface instead, it is possible to have an interface 
which without disorder moves for arbitrarily small driv- 
ing fields. The dynamics of the interface is then well 
defined for any strength of the disorder. However the 
most important observation is that now the conventional 
T = dynamics in which a spin flips if this lowers the 
energy of the system can be generalized to finite temper- 
atures. Since for weak disorder only small driving fields 
are needed at low temperatures conventional Glauber dy- 
namics can be used without running into the problem 
of spontaneous domain growth in the metastable phase. 
For small, bounded disorder there occurs a separation of 
time scales in the sense that for the system considered 
within the necessary simulation time only interface mo- 
tion is observed. The time scale for which the system 
runs into thermal equilibrium - which Glauber dynamics 
necessarily does of course - is many orders of magnitude 
larger, for more details see below. This clear separation 
of time scales makes a numerical study of temperature 
effects close to the pinning transition feasible. 

The random-field Ising model we consider is defined by 
the Hamiltonian 

n = - j J2 - H J2 a >-J2 H ^ w 

<ij> i i 

where Cj = ±1 are Ising spins on a two dimensional lat- 
tice. The first sum describes the ferromagnetic nearest 
neighbor interaction (J > 0). The random fields Hi are 
taken from a distribution which is constant within an 
interval [—A : A] and hence are bounded. H is the ho- 
mogenous driving field. 

The lattice we consider is shown in Fig. 0a). It has 
the structure of a square lattice but it is rotated by an 
angle of 7r/4. Hence, the boundaries and the initially 
flat interface of the system are in (1 Indirection. We use 
periodic boundary conditions parallel to the interface. 
Therefore all spins at the initial interface have exactly 
two of its four nearest-neighbor bonds broken, i. e. it 
does not cost any exchange energy to flip a spin at the 
interface. Hence, without disorder there is no pinning of 
the interface. 




FIG. 1. a) (upper picture) Schematic picture of a domain 
wall in the simulated lattice configuration, b) (lower picture) 
Domain wall of a simulated system of size 100 x 100 at the de- 
pinning field (H c = 1.29J) for finite temperature (T = 0.1J). 

In the direction perpendicular to the interface we use 
fixed boundary conditions. Spins in the lowest line are 
always up (in direction of the driving field) and in the last 
line are always down. At the beginning of the simulation 
all spins except of the lowest line are down and, hence, 
the interface is above the first line. Later, during the 
simulation the interface moves up due to the driving field. 

For finite temperature we perform Monte Carlo simu- 
lations where the spins are updated randomly using a 
heat-bath algorithm. For zero temperature this algo- 
rithm naturally crosses over to an algorithm where spins 
are only flipped when its energy decreases when flipped. 

Fig. [j]b) shows a system of size 100 x 100 spins dur- 
ing the simulation. The strength of the random field 
is A = 1.5 J and we will use this value throughout the 
whole paper. The driving field is H — 1.29J. We will 
show later that this value is the critical depinning field 
for A = 1.5 J at zero temperature. The interface seen is 
not pinned due to a finite temperature of T = 0.1J used 
in this simulation. No thermal fluctuation or domains 
appear in the unstable phase (upper part). This is im- 
portant for the following reasons. In equilibrium the two 
dimensional random-field system has no long range or- 
dered phase JTg) ] . Since Monte Carlo simulation in princi- 
pal leads to equilibrium properties one could expect that 
the system splits into domains spontaneously so that the 
concept of a single domain wall within the system is no 
longer useful. But this is not the case due to a sepa- 
ration of time scales mentioned above. Within the sin- 
gle spin flip method the growth of a domain must start 
with flipping one spin. The minimum energy needed for 
this process is AE = 2(4J — H — A) which is 3.42J 
here for the critical driving field. The corresponding 
flipping probability within a Monte Carlo simulation is 
exp(-A£/T) w 7 x 10" 14 for T = 0.1J. Therefore, 
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a thermal fluctuation within the bulk is far from being 
possible within our simulation time which is less than 10 5 
Monte Carlo steps per spin (MCS) - it happens on much 
larger time scales. Hence, as long as we restrict ourselves 
to low enough temperatures we can perform simulations 
without the possibility of spontaneous domain growth in 
the metastable phase. 

The velocity of the domain wall we define as the time 
derivative of the magnetization in a steady state. Note 
that at the beginning of the simulation the wall is flat 
and consequently the system is not in a steady state. In 
the steady state the magnetization grows linearly with 
time with fluctuations from sample to sample, of course. 
The velocity is obtained as averaged slope in this linear 
region. We average the velocity over many systems (10- 
160, depending on system size and how close the value of 
the driving field is to the critical one) and perform simu- 
lations for different system sizes (24 x 400 . . . 400 x 400) 
in order to investigate finite size effects systematically. 




FIG. 2. v versus H for T — and for finite temperatures 
ranging from 0.01J to 0.16 J. The system size is 400 x 400. 

Our results for the domain wall velocity v versus driv- 
ing field are shown in Fig. for different finite temper- 
atures and T = 0. During this simulation we consider 
for each temperature only driving fields which are large 
enough so that the wall crosses the system within roughly 
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FIG. 3. v versus H — H c with H c = 1.29 J determined from 
finite size scaling (see Fig. 4). T — 0. 



2000MCS. Otherwise we interrupt the simulation. As 
Fig. [^suggests we get a nice depinning transition at H c w 
1.29 J which is smeared out for finite temperatures. 

The wall velocity v depends on system size as shown 
in Fig. U for T — 0. Finite size effects are common at 
equilibrium phase transitions and there it is known that 
the most reliable values for the critical exponents are 
deduced from a finite size scaling analysis. We found 
that finite size scaling also works in the present case. For 
zero temperature the wall can be pinned and within these 
simulations we set the velocity of a single system to zero 
if the wall did not cross the system within roughly 20000 
MCS and we interrupt the simulation when more then 
one half of the simulated systems have velocity zero. In 
contrast to the situation at ordinary phase transitions 
here the effect of finite size is reversed: the smaller the 
systems are the stronger vanishes the order parameter. 
We use the finite-size scaling ansatz 



v(L,H) = L 
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where for the scaling function v we demand v(x) 



for x 3> 1 and — in contrast to ordinary phase transitions 
— v(x) = for x — > 0. 
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FIG. 4. Finite size scaling plot of Fig. ^. 

Fig. [| demonstrates that the finite size scaling ansatz 
above works well. Note that the scaling function has 
the unusual property mentioned above, i. e. it goes to 
zero for H — > H c . From the finite size scaling we get 
H c = 1.290J ± 0.003J, (3 = 0.35 ± 0.04, and v = 1.0 ± 
0.05. Alternatively, these results for /3 and H c can also 
be obtained directly from Fig. || by fitting those data 
which do not show finite size effects to a power law v ~ 
(H — Hey ■ The corresponding line is also shown in 
Fig. ^. Note also, that this value for j3 is in agreement 
with the earlier result in |Tl] | within the error bars. 

In order to analyse the finite temperature effects we 
first note that for the largest system size we investigated 
(400 x 400) there are no finite size effects for velocities v 
well above 0.2 as suggested by Fig. [| For L = 400 all 
data points except the one close to 0.2 are on a straight 
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line, i.e. increase as a power law without L-corrections. 
Hence, it is safe to neglect size effects for the finite tem- 
perature data in Fig. |[ For the smearing of the transition 
by temperature we expect a scaling behavior Q 

v(H, T) = T x ' 5 v ({H - HJT- 1 /^ . (3) 

A corresponding scaling plot of our data is shown in 
Fig. |. 
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FIG. 5. Scaling plot corresponding to Eq. ^ for the data of 
Fig. | 

We obtain a perfect data collapse with the following 
parameters: H c = 1.290J ± 0.005J, (3 = 0.33 ± 0.02, 
and 8 = 5.0 ± 0.3. To the best of our knowledge this 
is the first work in which the exponent 8 is determined 
for the system considered. Note that the perfect data 
collapse observed leads to exponents with only very small 
statistical errors. 

Finite size scaling and finite temperature scaling both 
give independent values for the critical field H c and for 
the exponent (3, respectively. It is rather satisfying that 
in both cases the numerical values agree within the error 
bars. 

As was mentioned in the introduction there are argu- 
ments in favor of the conjecture that the motion of a 
domain wall in a random-field system can be described 
by the EW equation with quenched disorder at the 
depinning transition. For this equation the scaling rela- 
tion u(2 — () — 1 has been derived jl2| by a functional 
renormalization group scheme in 4 — e-dimensions. This 
scaling relation has been claimed to be valid to all or- 
ders of e [jl3| . If we adopt this view we are therefore able 
to determine the roughness exponent of the random-field 
system without investigating the morphology of the in- 
terface. The value we obtain, C ~ F agrees with the value 
obtained from an e-expansion which is also claimed to be 
exact in all orders of e j[^,[L3) . Values of this exponent 
obtained numerically by integrating the EW equation or 
an automaton version of it |L4|] scatter between 0.7 and 
1.25 Jl^,|lJ]. From the scaling relation v{z -()=/? we 
get z sa 4/3 for the dynamic exponent which - as well 



as our result for (3 - is also in agreement with the re- 
sults of Nattermann et al. Interestingly, our value 
for the temperature exponent 8 fulfills the scaling rela- 
tion 8 = 2 + 1/(3 which has been derived by Tang and 
Stepanow p6[ by an extension of the functional renor- 
malization group scheme mentioned above to finite tem- 
peratures. All our findings support that the motion of a 
domain wall in a random-field system can be described 
by the EW equation. However, for an interface moving 
with a finite velocity also KPZ-like non-linearities could 
become relevant but from our numerical data we can- 
not extract any conclusions concerning a corresponding 
crossover of the values of exponents or scaling laws. 

To conclude, we investigated the influence of finite tem- 
peratures on the depinning transition of a driven [11]- 
interface by a Monte Carlo simulation of the two dimen- 
sional random-field Ising model. We used bounded disor- 
der and low temperatures in which case a clear separation 
of time scales occurs in the sense that no spontaneous 
growth of domains in the metastable phase appears. We 
derived the order parameter exponent j3 « 1/3 as well 
as the exponent v « 1 of the correlation length using 
finite size scaling. The corresponding dynamic exponent 
z ~ 4/3 and roughness exponent ( w 1 are determined 
via scaling relations. The exponent 8 describing the influ- 
ence of the temperature on the depinning of the domain 
wall which was unknown before was determined to be 
Spss 5. 
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